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A stability analysis of a spherically symmetric star in scalar-tensor theories of gravity is given 
in terms of the frequencies of quasi-normal modes. The scalar-tensor theories have a scalar field 
which is related to gravitation. There is an arbitrary function, the so-called coupling function, which 
determines the strength of the coupling between the gravitational scalar field and matter. Instability 
is induced by the scalar field for some ranges of the value of the first derivative of the coupling 
function. This instability leads to significant discrepancies with the results of binary-pulsar-timing 
experiments and hence, by the stability analysis, we can exclude the ranges of the first derivative 
of the coupling function in which the instability sets in. In this article, the constraint on the first 
derivative of the coupling function from the stability of relativistic stars is found. Analysis in terms 
of the quasi-normal mode frequencies accounts for the parameter dependence of the wave form of 
the scalar gravitational waves emitted from the Oppenheimer-Snyder collapse. The spontaneous 
scalarization is also discussed. 


I. INTRODUCTION 

General relativity is not the only viable theory of gravity. Scalar-tensor theories of gravity are among the alternative 
theories. The scalar-tensor theories contain a gravitational scalar field and an arbitrary function that determines the 
strength of the coupling between the scalar field and matter. With some choice of the coupling function, the scalar- 
tensor theories can pass all the present experimental tests of gravitation. From a theoretical point of view, it has 
been pointed out that the scalar-tensor theories arise naturally as the low-energy limit of string theory ||^,||. Until 
now almost all tests of high precision have been performed only in a weak field. It has been pointed out, however, 
by Damour and Esposito-Farese |^,|j that there are considerable deviations of predictions in a strong field given by 
some viable scalar-tensor theories from those given by general relativity. 

In the scalar-tensor theories of gravity, as in the case of general relativity, generation of gravitational waves is a 
typical phenomenon in a strong field and hence it is expected |^,|| that the projects (LIGO 0, VIRGO 1^, GEO j|] 
and TAMA [|^) of gravitational wave observations by laser interferometers will test the viable gravitational theories 
with high accuracy. Unlike general relativity, there is a scalar mode in gravitational waves in the scalar-tensor theories. 
Since this scalar gravitational wave can be detected by laser interferometers or resonant mass antennas and can 
be separated from tensor modes, it has crucial importance with respect to the observational tests of the theories of 
gravity. The scalar gravitational waves from the Oppenheimer-Snyder collapse, i.e., collapse of a spherically symmetric 
and homogeneous dust ball, have been calculated numerically as reported in Ref. |^, where behavior of the scalar 
field within the dust, the observed wave form of the scalar gravitational waves, and its dependence on the value of the 
first derivative of the coupling function are examined. From the observed wave form of the scalar gravitational waves 
we can obtain information concerning the first derivative of the coupling function. The dependence of the observed 
wave form on the first derivative of the coupling function can be explained by stability analysis of the quasi-normal 
modes since the behavior of the scalar field under the outgoing wave condition before the formation of the horizon is 
characterized by the quasi-normal modes that satisfy the regularity condition at the center. 

The stability analysis is motivated by the gravitational test through pulsar-timing experiments. By measuring the 
change of the period of pulsation we can obtain information regarding scalar “form factors” of the pulsar, and those 
form factors are related to the variation of the moment of inertia under the influence of the companion star Q]. It 
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was shown that in some class of the scalar-tensor theories there is a solution with non-perturbative strong-held effects 
and it has form factors that are signihcantly different from those in general relativity In order to obtain the 

correct gravitational theory, it is important to investigate the internal structure of highly relativistic stars and check 
the consistency with the results observational tests. The stability analysis reveals whether the deviation of the stable 
equilibrium solution of the relativistic star from that in general relativity remains small or not. If, with some choice 
of the coupling function, the stable equilibrium solution suffers from non-perturbative effects by the coupling between 
the scalar held and matter and has signihcantly different scalar form factors from those in general relativity, the 
results of the pulsar-timing experiments will exclude that choice of the coupling function. 

The condition for the stability of relativistic stars of a perfect huid in the scalar-tensor theories for radial and 
time-symmetric perturbations was examined by Bruckman and Velazquez |p^ . They found that the stability depends 
on the choice of the coupling function. 

This paper is organized as follows. In Sec. || we review the scalar-tensor theories of gravity and derive the basic 
equation. In Sec. Ill we introduce the quasi-normal modes of the scalar held and discuss some useful conditions 
concerning them. In Sec. IV we introduce stellar models. In Sec. ^ we give numerical results on stability versus 


instability for these stellar models. In Sec. 
units in which c = 1. We follow the MTW 
tensor. 


1 we summarize the results and discuss their implications. We use the 
]I3| sign conventions for the metric tensor, Riemann tensor and Einstein 


II. BASIC EQUATION 


A. Scalar-Tensor Theories of Gravity 


We consider a class of scalar-tensor theories of gravity in which gravitation is mediated by a long-range scalar held 
in addition to space-time curvature. The action in the Einstein frame formulation is given by |I4|] 


I = 


1 


IGttG 




( 2 . 1 ) 


where the metric tensor in the “unphysical” Einstein frame is related to the metric tensor in the “physical” 
Brans-Dicke frame by the conformal transformation as 


= A'^{(p)g 




( 2 . 2 ) 


Here, G* is the “bare” gravitational constant and Lp is the gravitational scalar held, i?* is the scalar curvature of 
and 'I'm represents matter helds collectively. The arbitrary function A{ip) is related to a coupling function which will 
be introduced below. The Brans-Dicke frame metric g^^i, is universally coupled to matter in the model we consider 
here, and hence the Einstein equivalence principle is valid. 

The held equations are given by 


□ = —47rG*a(i^)T*, 


(2.3) 

(2.4) 


and the equations of motion are given by 




(2.5) 


where 


T* 

a{ip) 


2 A'^{(p)g^:fj_,^] _ 

9 ^ ^ 9 * 111 ^ 

T ^ 

dhiA{‘^) 

d^p 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


In the above equations, -n, and are the Einstein tensor, covariant derivative and d’Alembertian of 

respectively, and is the stress-energy tensor of matter in the Brans-Dicke frame. 
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The function a(t p) is a coupling function between the scalar field and trace of the stress-energy tensor of matter, 
as seen in Eq. (|2.4|) . If a{}p) is constant, the theory reduces to the Brans-Dicke theory If a{ip) = 0 , the theory 
reduces to general relativity. Hereafter we consider an asymptotically flat space-time and we assume ^ Lpo at 
spatial infinity. We use units in which G* = 1, and we fix the freedom of the constant conformal transformation by 
requiring = 1. 

Here we briefly summarize the present constraints on the coupling function a{Lp) through the results of solar-system 
test experiments. The solar system is a laboratory in a weak gravitational field, and hence the parametrized-post- 
Newtonian (PPN) formalism works well. In the scalar-tensor theories introduced above, the so-called Eddington 
parameters among the PPN parameters are expressed by the coupling function a((/?) as [p^ 


1 — lEdd = 
Psdd — 1 = 


2ag 

l + «o’ 


/3oao 

2(l + ag)2’ 


(2.9) 

( 2 . 10 ) 


where 


ao = a(</5o), (2.11) 

Po = ^{^o), (2.12) 

and the values of the other PPN parameters that enter in the first post-Newtonian approximation are identical to 

those for general relativity . The results of the experiments of the light deflection constrain ■jEdd as 0 

bBdd = 0.9996 ±0.0017, (2.13) 


and this constraint reduces to 


al < 0.001. 

(2.14) 

The results of the lunar laser-ranging experiments constrain ^Edd as |]l8|] 


= 0.9998 ±0.0006, 

(2.15) 

and this constant reduces to 


Poal = -0.0004 ±0.0012. 

(2.16) 


B. Perturbation Equations 


From this point we restrict our attention to the scalar-tensor theories in which 


ao = 0. 


(2.17) 


This choice is completely consistent with the experimental constraints Eqs. ( 2.14 ) and dm)- The scalar-tensor 
theories in which ao = 0 pass all weak-field tests of solar-system experiments that have been carried out until now, 
for an arbitrary value of Pg. Further by this assumption we can seek the pure effect of the first derivative /3o of the 
coupling function a{ip) on the sta bilit y of a star . 

For ao = 0, the field equations (2.3) and (2.4) allow the solution 


a - 
'E — 71(B) 

^ = V^O, 


(2.18) 

(2.19) 

( 2 . 20 ) 


where the set of g\iv and Tpii, is a solution of the Einstein equation, i.e.. 




( 2 . 21 ) 
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and therefore 


( 2 . 22 ) 


^(E)TiE> = 0. 


Here, and are the covariant derivative and Einstein tensor of respectively. In this solution the Einstein 
frame expression agrees with the Brans-Dicke frame expression: 


T — T 

fiiy — • 

If we consider linear perturbations of this solution, the coupled differential equations for the perturbations decouple 
to an equation for the scalar field perturbation and equations for the metric and matter perturbations. This is because, 
for ao = 0, the perturbation of the conformal factor A‘^{(p) defined in Eq. (2.2) by the scalar field perturbation from 
the constant value (po is a second-order infinitesimal. For the same reason, the perturbations and of g^^i, 

and in the Einstein frame are the same as and of g^i, and in the Brans-Dicke frame up to linear 
order: 


= 9\S) + ^9iiv + higher order terms, 

^ + higher order terms, 

ip = ipo + 6ip. 


From Eqs. (^), (^) and (|2.5|), we obtain the field equations for the linear perturbations as 


and the equations of motion as 


SGf,, = SttST^,, 


= 0 , 


(2.23) 

(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 


where is the d’Alembertian of gl3^■ Equations ( 2.26| ) for the metric and matter perturbations are identical to 
those for general relativity. Since we are now interested in instability induced by the gravit ationa l scalar field we do 
not consider these equations here. Thus it is sufficient for our purposes to consider only Eq. ( 2.27 ) for the scalar field 
perturbation. In what follows we omit the symbol (E) to simplify the notation. 

We assume that the unperturbed solution is static and spherically symmetric. Then, without loss of generality, the 
metric is written in the simple form 


ds^ = -k -h r‘^{de‘^ + sin^ Odcf^). 


(2.29) 


In this space-time, Eq. ( |2.27|) becomes 


_2<i,d‘^5ip , 

Q 



dt^ ' 


dr ) 

1 

■ 1 d 


1 d'^dip 

~ j,2 

sin 0 do 

86 j 

' sin^ 9 8(jP _ 


- 4:TrPoT6p> = 0 . 


We concentrate on perturbations with sinusoidal time dependence, and define ipujimir) as 




(2.30) 


(2.31) 


l,m 


where we have followed the sign convention of Ref. 0. Hereafter we abbreviate puiimir) as '^{r). Then from Eq. ( |2.30| ) 
we obtain a Schrodinger-like equation as 


d^p 

drl 


+ - E ( r *)]'0 = 0 , 


(2.32) 


where we have defined the “generalized tortoise coordinate” r* as 
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dr^ = '^dr. 


(2.33) 


The effective potential, V{r^,), is defined as 


V{r^) ^ - d^/JoTe^^, 


(2.34) 


where the prime denotes d/dr and Pq has been defined in Eq. ( ^.12| ). On the right-hand side of Eq. ( 2.34 ), the first 
term is a “curvature potential”, the second is a “centrifugal potential”, and the third comes from the fact that the 
coupling function a{ip) depends on the scalar field tp, and hence a{ip) is not constant when the perturbation of the 
scalar field (f is present. 


III. STABILITY ANALYSIS 


A. the Quasi-Normal Modes 


For any complex number tu, a general solution of Eq. (2.32) consists of a linear combination of two independent 
solutions as 


+ C (r*), 

where and 'ip~ are purely ingoing and outgoing waves. That is, 

in the limit r* —> oo. The coefficients C^(io) and C~{u) are determined by 


(3.1) 


(3.2) 

(3.3) 


C+{u,) 

c-{^) 


1T(V’+,^/’-)’ 

1T(V’,V’+) 

1T(V’-,V'+)’ 


where W (/, g) is the Wronskian defined as 


W{f,g) 


dg df 

dr^ dr^ 


(3.4) 

(3.5) 

(3.6) 


For a general complex number uj there may be no solution that satisfies the regularity condition {pj = 0) at the 
origin and outgoing wave boundary condition at infinity = 0) simultaneously. For some special, discrete 

values of the frequency oj there exists a purely outgoing wave solution that satisfies the regularity condition. We 
have special interest in this solution as an eigenmode of the scalar field. We call the solution a quasi-normal mode 
and this special complex number io a quasi-normal mode frequency. The physical meaning of the quasi-normal mode 
of the scalar field is that when the scalar field is perturbed by some effect other than the incidence of scalar waves, 
the scalar field oscillates with radiation of scalar waves, and its characteristic oscillation modes are described by the 
quasi-normal modes. 


B. Stability of the Quasi-Normal Modes 


With the sign convention in Eq. (2.31), if the imaginary part of the quasi-normal mode frequency uj is positive, the 
quasi-normal oscillation is damped, and its damping rate is given by this imaginary part. If the imaginary part of lo 
is negative, however, the quasi-normal mode is unstable, and the amplitude of the perturbation grows exponentially 
in the growth rate given by the absolute value of the imaginary part. The unstable quasi-normal mode develops from 
a regular initial set of data, because the amplitude of the perturbation decreases exponentially to zero in the limit 
r* —> oo. Whether the quasi-normal mode is stable or not is thus determined by whether its complex frequency is in 
the upper half plane or in the lower half plane. 
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As for the stability of the quasi-normal modes of Eq. ( 2.32| ), the following conditions, similar to those for the 
quasi-normal modes of spherical nonrotating stars in general relativity 1^ , hold: (i) An unstable quasi-normal mode 
of any I has a purely imaginary frequency, i.e., an unstable mode grows exponentially in time without oscillation, 
(ii) If the effective potential V is non-negative everywhere, there is no unstable quasi-normal mode, (iii) Consider 
a sequence of effective potentials smoothly parametrized by a real variable A. Then there are discrete quasi-normal 
mode frequencies, and each quasi-normal mode frequency, denoted by ujn{X), moves smoothly along a trajectory in 
the complex plane as A varies smoothly, (iv) A quasi-normal mode frequency is bounded, (v) Consider a sequence of 
effective potentials smoothly parametrized by a real variable A. If there is a critical value for any quasi-normal 
mode frequency LUn such that Im[a;„(A)] > 0 for A > Ajf** and Im[a;n(A)] < 0 for A < Ajf**, then a;„(Ajf**) = 0. In 
other words, when the stability of a quasi-normal mode changes as A varies smoothly, its frequency is zero. 

Condition (i) is the direct result of the fact that the finiteness of the norm of an unstable quasi-normal mode and 
hermitianity of the derivative operator 


drl 


■E, 


(3.7) 


guarantee that the eigenvalue is real for the unstable quasi-normal mode. Condition (ii) is equivalent to the 
absence of bound states in the non-negative potential in quantum mechanics. The proof of condition (iii) is almost 
the same as that for the quasi-normal mode frequencies of gravitational waves in general relativity The sketch 
of the proof is as follows. First we take u = in place of ^ (cf. Appendix E- The regularity condition 

requires du/dr = 0 at the origin, and we fix the normalization by requiring m = 1 at the origin. The solutions u, 
and therefore tp depend smoothly on A because each term in the differential equation ( |2.32[ ) depends smoothly on A. 
So A), which is defined by C+(a;) for A, is a smooth function of A. The solutions u and ther efore ip are also 

analytic functions of uj because the analyticity of each term for uj in the differential equation (2.32) guarantees the 
Cauchy-Riemann relations for its solutions. Therefore A) is an analytic function of uj. From the property of an 

analytic function, any zero of C'^{uj] A) for any given A, denoted by ujn[\), is isolated, and the number of the zeroes is, 
at most, countably infinite. If we were to assume the discontinuity of each quasi-normal mode frequency oJnW with 
respect to A, a contradiction would result with the smoothness of A) with respect to A through the maximum 

principle. Therefore we conclude that each quasi-normal mode frequency tCn(A) varies smoothly with respect to A. 
Condition (iv) can be easily shown by observing that in the limit |w| ^ oo, the solution regular at the origin is not a 
purely outgoing wave. Condition (v) is the combined result of conditions (i), (iii) and (iv). 

We define a critical value as a value of Po such that for Pq > (<)/3g’’*‘ there is no unstable mode and for 
Po < {>)Po^^* there is an unstable mode. The corollary of condition (ii) is that if T = < 0 in the stellar interior, 

which is considered to be satisfied in most of the physical situations, then there is a critical value and for any 

Po larger than this value there is no unstable quasi-normal mode. The corollary of condition (iii) is that each quasi¬ 
normal mode frequency moves smoothly in the complex plane as we change Po smoothly. The corollary of condition 
(v) is that there is a zero-frequency quasi-normal mode at Po = /Jg’"**. 


IV. MODELS 


To obtain further understanding of the instability induced by the scalar field we consider stellar models that describe 
relativistic stars. First we take as a stellar model the equilibrium solution of an incompressible fluid (see Box 23.2 of 
MTW |13). In this case the matter is a perfect fluid: 


= pu^u^ + + Ufj,Uu). 

The trace of the stress-energy tensor is given by 

T = Ti^^ = -p + 3P. 

The energy density and pressure are given by 


p(r) = / ^0 (0 < r < R) 


P{r) = 


_ ) Po 


0 (i? < r) 

3(i_^)V2_(i_^)1/- 


(0 < r < i?) 
{R < r) 


(4.1) 

(4.2) 

(4.3) 

(4.4) 













where M and R are the gravitational mass and radius of the star, 


M = 



The metric is given in the form of Eq. ( 2.29 ), where 

3 [3(1-)'/"-(! 


1 _ 2M 


g2'I-(r) 

m{r) 


2m(r) 


( ^irpor^ {0 < r < R) 
[M (i? < r) 



{0 < r < R) 
{R < r) 


(4.5) 


(4.6) 

(4.7) 

(4.8) 


There exists a smallest radius of a star of homogeneous density. From this, we have the condition 

R 9 
M ^ 4’ 


(4.9) 


We can easily see that the negativity of T in the stellar interior holds if and only if 


i? m 

and the dominant energy condition holds if and only if 

R 8 
M ^ 3' 


(4.10) 


(4.11) 


Next we introduce a pol ytrop ic stellar model. The metric for a static and spherically s ymm etric space-time is 
written in the form of Eq. ( |2.29|) . The stress-energy tensor of a perfect fluid is given by Eq. (|4.l| ). We integrate the 
Tolman-Oppenheimer-Volkoff equation 


dP .m + Airr^P 

dr ^ r[r — 2m] ’ 


where m{r) is defined by Eq. (^^) and determined by 

m(r) = 47r / p{r')r'^dr'. 

Jo 

The function <i>(r) is determined by integrating the differential equation 


dr 


1 dP 
P + p dr ^ 


(4.12) 


(4.13) 


(4.14) 


and the requirement of matching to the Schwarzschild space-time at the stellar surface. We consider the following 
polytropic equation of state: 


p = nrub + 


KriQinb ( n 


r - 1 V no 

r 
n 


P = Knorrib — 
\no 

rub = 1.66 X 10“^^g, 
no = O.lfm”^. 


(4.15) 

(4.16) 

(4.17) 

(4.18) 


We then take the parameters T = 2.34 and K = 0.0195 Q to fit the more realistic equation of state for high-density 
nuclear matter 
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FIG. 1. Z-dependence on the effective potential. The ordinate is the effective potential and the abscissa is the areal coordinate. 
The stellar radius R is fixed to lOM. The theoretical parameter f3o is fixed to 0. 



r/M 

FIG. 2. /3o-dependence on the effective potential. The ordinate and abscissa are the same as in Fig. 1. The stellar radius R 
is fixed to lOM. The angular momentum I is fixed to 0. 
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FIG. 3. Comparison between the effective potentials of an Z = 0 mode for /3o = —5 of the incompressible fluid model and 
polytropic model. The ordinate and abscissa are the same as in Fig. 1. The solid and dashed lines denote the incompressible 
fluid model and polytropic model respectively. In the polytropic model, the total central density is 1.16 x 10^®g/cm^, and the 
stellar mass is l.TdM©. The stellar radius is 4.94M in both models. The details of the polytropic stellar model are given in 
the Sec. 0 and Appendix]^ 





































The effective potentials of the stellar model of an incompressible fluid are shown in Figs. 1, 2 and 3. Figure 1 
shows the Z-dependence of the shape of the effective potential of the incompressible fluid model for R = lOM. The 
parameter /3o is fixed to 0. The centrifugal potential raises signihcantly the total effective potential for I > 1 modes. 
Therefore, from condition (i), an Z = 0 (spherically symmetric) mode is more likely to be unstable than any Z > 1 
mode, and therefore we concentrate our attention to the Z = 0 mode. Figure 2 shows the /3o-dependence of the shape 
of the effective potential of the incompressible fluid model for Z = 0 and R = lOM. The term which comes from the 
coupling between the scalar field and matter depends on /3o. As /3o becomes smaller, the quasi-normal modes become 
more likely to be unstable because the potential well becomes deeper. In particular, from condition (i) and Fig. 2, 
no unstable quasi-normal mode exists in the theory in which /3o = 1 for i? = lOM in this model. Figure 3 shows the 
effective potential of the stellar model of an incompressible fluid and that of the polytropic stellar model of an Z = 0 
mode for Pq = —5. In Fig. 3 we hx the central total energy density to be 1.16 x lO^^g/cm^ in the polytropic stellar 
model and obtain the stellar radius 12.7km, mass 1.74M0 and radius-to-mass ratio 4.94. We plot, for comparison, 
the effective potential of the model of an incompressible fluid whose radius-to-mass ratio is also fixed to 4.94. The 
details of the numerical calculation of the polytropic stellar model are described in Appendix 


V. NUMERICAL RESULTS 


To search for we integrate the differential equation (2.32|) for w = 0 from the origin under the regularity 


condition. The regularity condition is given by the Taylor expansion of the solution around the ori gin. The expansion 
of the regular solution around the origin is described in Appendix After we integrate Eq. ( ^.32[) to the stellar 
surface we calculate the Wronskian with an exterior solution that is regular at r = oo. As shown in Appendix]^ we 
know that the exterior solution for to = 0 that is regular at r = oo is given by the hypergeometric function. Then 
we examine whether or not the Wronskian vanishes. However, the existence of a zero frequency quasi-normal mode 
does not necessarily imply the onset of instability because condition (v) does not deny the existence of a for 

LOn such that = 0 but Im[w„(A)] > 0 for A > and also for A < In order to check the onset of 

instability we relax the assumption w = 0. We consider a pure imaginary frequency, i.e.. 


OJ = —ic. 


(OO) 


(5.1) 


and under the regularity condition we integrate Eq. ( ^.32 ) from the origin to the surface and thereafter from the 
surface to the asymptotic region \uj\r ^ 1 and r M. Then we calculate the Wronskian with the asymptotic 
expansion of ip~ and the coefficient The details of the asymptotic expansion are described in Appendix 

When Po is changed to a slightly smaller (larger) value than the value for the emergence of a zero-frequency mode, 
we can confirm that an unstable quasi-normal mode emerges by observing that C~^{u! = —ic) crosses zero as c is 
increased from zero. 



FIG. 4. Stability versus instability. The ordinate is the theoretical parameter /3o and the abscissa is the stellar radius. The 
solid lines denote the lines Po = Po"''*', and the dashed lines denote the onset of instability of another quasi-normal mode in 
the incompressible fluid model. The open circles denote the points Po = Ptf''*', and the filled circles the onset of instability of 
another quasi-normal mode in the polytropic model. 


9 




















60 


02 . 


40 

20 

0 

-20 

-40 

-60 


-unstable 


stable 


unstable 

• • 


3 3.5 4 4.5 5 5.5 6 6.5 7 

R/M 


FIG. 5. Magnification of the region of smaller stellar radius in Fig. 4. 


TABLE I. The critical values for the polytropic stellar models, pc, M and R are the total energy density at the center, 
gravtitational mass and radius, respectively. Only the equilibrium solutions that satisfy the necessary condition for stability 
dM/dpc > 0 are in the table. 



Pc{g/cm^) 

M/Mq 

i?(km) 

pcr^t 

2.26E+15 

1.95 

11.2 

-4.90 

1.92E+15 

1.93 

11.6 

-4.59 

1.45E+15 

1.85 

12.2 

-4.35 

1.16E+15 

1.74 

12.7 

-4.39 

1.03E+15 

1.67 

12.9 

-4.48 

8.96E+14 

1.57 

13.2 

-4.67 

7.70E+14 

1.45 

13.4 

-4.97 

5.32E+14 

1.12 

13.7 

-6.21 

3.10E+14 

0.670 

13.7 

-9.94 

2.04E+14 

0.410 

13.3 

-15.5 

l.OOE+14 

0.161 

12.3 

-36.2 
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The results of numerical calculations are shown in Figs. 4, 5 and 6. Figures 4 and 5 are diagrams of stability 
versus instability in the Po vs. R/M plane. The solid and dashed lines denote the emergence of the first and second 
unstable modes for the incompressible fluid model. Here we have restricted our attention to the stellar model satisfying 
R > (8/3)M so that the dominant energy condition may be satisfied. If we consider more realistic models, there is the 
necessary condition, dM/dpc > 0, for stability, where pc is the total density at the center. From this condition, the 
ratio R/M cannot be smaller than some critical value that depends on the equation of state, for example, 3.18 — 4.63 
for models in Ref. [^. If the radius is smaller than this critical radius, the stellar model is already unstable even if 
Po = 0, i.e., even if there is no coupling between the scalar field and matter. The unstable region above the line /?o = 0, 
as seen in Figs. 4 and 5, is due to the violation of negativity of T = —p + 3P for the stellar radius R < (18/5)M. As 
seen in Fig. 4, for R'^ M the curve of /3o = P/p''* can be fitted by the straight line 


ocrit 

Po 


TT^ R 
12 M’ 


(5.2) 


and this line is derived by the emergence of the unstable mode in the square well potential ignoring the space-time 
curvature and matter pressure, and assuming a homogeneous density distribution. The open and filled circles denote 
the emergence of the first and second unstable quasi-normal modes for the polytropic stellar models that satisfy the 
necessary condition for stability, dM/dpc > 0. They are also shown in Table I. Figures 4 and 5 suggest that the 
critical value /Jg’"** for the emergence of the first unstable mode has little dependence on the equation of state, as a 
function of the ratio R/M. Figure 6 displays the emergence of an unstable quasi-normal mode for the incompressible 
fluid model of i? = lOM. Since the eigenvalue of the first unstable mo de is approximately given by the depth of the 
effective potential V and the third term dominates other terms in Eq. ( ^.34 ), if Pq is considerably smaller than P^ 
the growth time r for the first unstable mode is approximately given by 


crit 
'0 ’ 


r ~ 




(5.3) 


where tq denotes the free-fall time of the star. 
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FIG. 6. Emergence of an unstable quasi-normal mode. The ordinate is = —ic) and the abscissa is c = —Im[a;]. 



VI. SUMMARY AND DISCUSSIONS 

In the scalar-tensor theories of gravity in which ao = a:((po) = 0, the field equations allow a solution that consists 
of a constant scalar field and a solution of the Einstein equation for metric tensor and matter. The field equations for 
linear perturbations of this solution are decoupled to an equation for the scalar field perturbation and equations for 
the metric and matter perturbations. The equations for the metric and matter perturbations are identical to those 
in general relativity. 

The stability analysis of the equation for the scalar field perturbation is given in terms of the quasi-normal mode 
frequencies. From a general argument on the quasi-normal modes of the scalar field in the static and spherically 
symmetric space-time, the following results turn out to be true: An unstable mode grows monotonically. The presence 
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of a non-negative effective potential implies the absence of unstable quasi-normal modes. When a quasi-normal mode 
changes its stability, its frequency is zero. 

From the shape of the effective potential, we find the following tendencies with regard to stability: As I becomes 
smaller, the quasi-normal modes of the scalar field become more likely to be unstable . Stability of the quasi-normal 
modes also depends on the first derivative /3o of the coupling function a{(p) between the scalar field and matter. As 
Po becomes smaller, the quasi-normal modes of the scalar field become more likely to be unstable if T = —p + 3P < 0 
in the stellar interior. 

In order to seek for the critical value for the onset of instability of a quasi-normal mode of the scalar field, we 
examined whether or not a zero-frequency exists, because the frequency of the marginally stable quasi-normal mode 
that changes its stability must be zero. Then we confirmed the onset of instability by observing the emergence of an 
unstable mode. 

Using the incompressible fluid model and polytropic model of a star we found the critical value numerically 
as a function oi R/M. If i?/M ~ 4 for neutron stars and if these stellar models approximate such neutron stars 
well, our numerical results suggest that a constraint on the parameter /3o would be /3o ^ —5 from the stability of the 
neutron stars. In the theories in which oq = 0 and Po ^ ~5 (this value depends slightly on the equation of state for 
high-density nuclear matter), the coupling between the scalar field and matter may change the stability of a stellar 
solution that consists of a solution for the Einstein equation and a constant scalar field. Even if a stable equilibrium 
solution exists for Po < —5, such a solution must contain non-perturbative effects due to the coupling between the 
scalar field and matter and it will have significantly different form factors for the pulsar-timing experiments from 
those in general relativity Q. Even for non-relativistic stars, such as white dwarfs, main-sequence stars, planets and 
so on, there is a critical value /3g''*‘ ~ —R/M for their stability. 

Damour and Esposito-Farese discovered that, for a coupling function of the form 

a{ip) = Pp, (6.1) 


when P < —4, non-perturbative effects develop for massive neutron stars. That is, the star has a nontrivial finite 
“scalar charge” .and the scalar field displays nontrivial configuration even when the asymptotic value po is extremely 


small or zero The solution with the non-perturbative effects is significantly different from that in general 

relativity. These results raise a paradoxical problem in the limit tpo —>■ 0. If tpo = 0, both the field equations and 
boundary condition possess symmetry under the reflection transformation ip — ipo —> —((p — ipo), and hence they 
have the trivial solution (p{x) = ipo = 0, and the matter and metric tensor satisfy the Einstein equation in general 
relativity. It seems that this fact implies a discontinuity in the sequence of the solutions against ipo. Recently it has 
been discussed that the emergence of the non-perturbative effects on the equilibrium solutions of neutron stars is 
interpreted as a “spontaneous scalarization” analogous to a spontaneous magnetization of the ferromagnets [Q. The 
critical value ~ — 5 for the onset of instability of a quasi-normal mode for relativistic stars obtained here by our 
analysis is approximately equal to the critical value in which an equilibrium solution with non-perturbative effects 
for relativistic stars emerges for og = 0 In the context of spontaneous scalarization, we have shown that if 

Po > the equilibrium solution with symmetry for the reflection transformation p — ipQ —{ip — pp) is stable, 

but if Po < P/p''* the solution with the symmetry is no longer stable. This fact strongly suggests that a spontaneous 
symmetry breaking occurs at Pq = PP"^*. It is still important to examine the stability of the equilibrium solution with 
non-perturbative effects. 

It was shown in Ref. that in the Oppenheimer-Snyder collapse the scalar field within the dust decays with 
oscillation for non-negative values of Po (see Fig. 14 of Ref. 0 ) and grows without oscillation for a large negative 
value of Po (see Fig. 15 of Ref. 0 ). This leads to the /3o-dependence of the observed wave form of the scalar 
gravitational waves, as seen in Figs. 10, 11 and 12 of Ref. |^. The results of the stability analysis given above 
account for the behavior of the scalar field in terms of the quasi-normal mode frequencies. For Po > /Jg”** no unstable 
quasi-normal mode exists, and only the damped oscillation modes are allowed, but for Po < PP^^ unstable modes exist 
and grow exponentially without oscillation. Those purely outgoing wave modes characterize the time development of 
the scalar field within the dust and therefore the observed wave form of the scalar gravitational waves. For Pq > PP'^*, 
the quasi-normal mode oscillations appear in the observed wave form before the last quasi-normal mode ringing of 
the formed black hole, as seen in Fig. 11 of Ref. 0 |. For Pq < Pp''^, the growing quasi-normal mode appears in the 
wave form before the last quasi-normal mode ringing of the black hole, as seen in Fig. 12 of Ref. [pd]| . 

Although we have assumed implicitly that the unperturbed solution has a regular origin, we can extend the results 
obtained in Sec. Ill to a solution that has an event horizon. In this case we have only to replace the regularity condition 
at the origin r = 0 with the ingoing wave condition at the horizon, r* ^ —oo. If we consider a Schwarzschild or 
Reissner-Nordtstrom black hole, coupling between the scalar fileld a nd m atter does not exist because T = 0 for the 
vacuum and electromagnetic field. Therefore the effective potential (2.34) is identical to that for a non-gravitational 
scalar field in general relativity. There are, however, some exotic black hole solutions in which T ^ 0, and the stability 
of such black holes may be altered due to the coupling between the gravitational scalar field and matter. 
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Although we have assumed ao = 0, the stability analysis given here applies to the perturbation equations truncated 
up to zeroth order in ao expansions around general relativity |lll , if oq 7^ 0. When we need to consider the effect of 
the nonzero value of oq seriously, we must treat the coupled equations for the perturbations of the scalar field, metric 
tensor and matter fields. The assumption that ao is extremely small seems to be reasonable if the first derivative /3o 
of the coupling function a{(p) is positive or if the scalar field has a small mass rriip ( 10®km < h/rriip < ) and the 

location (pm of the minimum of the potential of tp , V(p) = (l/87r)m^((^ — PmY', coincides with a zero of the coupling 
function a{p)^ because in these cases the cosmological attraction mechanism is effective and ao is attracted toward 
zero 
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APPENDIX A: EXPANSION OF REGULAR SOLUTION AROUND THE ORIGIN 


To suppress numerical errors in the stellar interior we integrate not tp but u defined as 

ip = 


(Al) 


Then the differential equation for u is 


($' - 'k') + 


2{l + l) 




In the exterior vacuum region the differential equation for ip is 

2M\~^ 2M 


-(e^^ - 1) + 47r/3oTe^ 


It = 0. 


(A2) 


r + [i - 


+ 1- 


T) 




1 - 


2M' 

r J. 


V’ = 0. 


(A3) 


In order to obtain the expansion of the solution u around the origin, we expand T, <I> and 'k around the origin as 


T = To + Tar^ + H-, 

(A4) 

$ = $0 + $2r^ + H-, 

(A5) 

\k = ^'o + 'k2r^ + 'k4r"‘ H-. 

(A6) 

If we take the incompressible fluid model, the coefficients are given by 


T 93^-2 

(A7) 

*’'(31-1)='’"’ 

(A8) 


(A9) 

47r 

^ “ 3(37-1)^°’ 

(AlO) 

8 2 37 - 2 2 

*‘=9' (37-I)"""' 

(All) 

0 

II 

. ^ 

(A12) 

4 

^2 = gTrpo, 

(A13) 

T 16 2 2 

^'4 = Po. 

(A14) 
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where we have defined 


7 = 



and 1/3 < 7 < 1, as seen from Eq. (4.9). 

Using these expansions, the regular solution u is expanded around the origin as 


00 

u = ^ ayr'", 

i /=0 


where the coefficients {oj,} are determined using the coefficients in Eqs. (R3), (A5) and (A 6 ) as 


oi = 03 = • • • = 0 , 

02 = + 21^2 - 21(1 + 2)^2 + iTTpoTo] oq , 

04 = (^-2 - $ 2 ) + 4Z($4 - 4 - 4 ) 

-2i{i + i)(4'4 + 4 - 2 ) + 4^/3o(r2 + 2ro^2)] 00 
+ + 2(/ + 2)^2 — 2(Z^ 21 2)^2 “t“ 47r/3o^o] ^ 2 } 

= 8 ( 2 i + 3 K 2 i + 5 ) f + (^ + 2 ')*. - + 21 + 2 )*, 

+AttPoTo] + 2Z$2 - 2l{l + 2)^-2 + AirpoTo] 

-2(21 + 3) [2L0^e-^'^°i^2 - $ 2 ) + 4Z($4 - ^' 4 ) 

- 21(1 + l)(«-4 + 4 - 2 ) + 4^/3o(r2 + 2ro^2)] } oo, 

0-6 ‘ ‘ ‘ ; 


and so on. 


(A15) 


(A16) 


(A17) 

(A18) 


(A19) 


(A20) 


APPENDIX B: NUMERICAL CALCULATION OF THE POLYTROPIC STELLAR MODEL 


We make the variables dimensionless: 


~ _ n 


P = 

P = 

f = 
rh = 


no O.lfm 
P 


-3 > 


nomb 

1.66 X 10i4gcm-3’ 

P 

P 

nomb 

1.49 X 10^®gcm“^s“2 

r 

r 

(nomb)~ 

■ 1/2 9.02 x106cm’ 

m 

m 


(noTOf,)“^/^ 1.22 X lO^^g 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


Hereafter we om it tilde s for s impli city. Then we obtain the following ordinary differential equations for n , m and 4* 
from Eqs. (HI), (PI) and (PI ): 


dn 

dr 

dm 

dr 

d^ 

dr 


n(P + p) m + Anr^P 
TP r(r — 2m) ’ 

dTTpr^, 

m + Anr^P 
r(r — 2m) ’ 


(B 6 ) 

(B7) 

(B 8 ) 
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where 


P 

p = n+ — 
P = Kn^. 


1 ’ 


We integrate Eqs. (|B6|), (B7) and (|BS|) 


from r = 0. The initial values of n, m and are 


n = He, 
m = 0, 
$ = 0 . 


(B9) 

(BIO) 


(Bll) 

(B12) 

(B13) 


In order to match the internal solution to the exterior Schwarzschild space-time at the stellar surface, we modify the 
calculated value $caJc to the true value of ^tme by 




$caZc(r) - ^iR)calc + 2 



(B14) 


To guarantee regularity at the origin, we use the following expansions of the solution around the origin: 


n = nc + n 2 r + n^r + ■ ■ 
m = + • • • j 

$ = $2r^ -I- -I-, 

p = Pc + P2r^ H-, 

p = p^+p^r'^ + py + ■ 

Then the coefficients are given by the following relations. 

Pc = Knl, 


Pc = nc + 


Pc 


r- 1’ 


W3 = gTrpc, 

n2 = + 47rPc), 

n2 

P 2 = TPe — , 

Uc 


P2 =n2 + 


P 2 


r-V 


W5 = ■r7rp2, 
5 


$2 = — 


P 2 


Pc + Pi 
1 


714 = 


4rp, 


772 (Pc + Pc) + 77c ( P2 - ^P2 ) f- (7773 + 47rPc) 


-|-77c(Pc -I- Pc)( 7775 -I- 47rP2) + 277c(Pc “t" Pc){m3 + 47rPc)7773] 
2 

+ rPe 


P 4 = -r(r-i)p, 


<1)4 = - 


P 4 




77c 

^ 2(^2 + P 2 ) 


Pc + Pc ‘^{Pc + Pc)^ 


(B15) 

(B16) 

(B17) 

(B18) 

(B19) 


(B20) 

(B21) 

(B22) 

(B23) 

(B24) 

(B25) 

(B26) 

(B27) 


(B28) 

(B29) 

(B30) 


APPENDIX C: STATIC SOLUTION IN THE EXTERIOR SCHWARZSCHILD SOLUTION 

If we assume the scalar field ^ does not depend on t, Eq. (^.32|) reduces to 
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Cp-lp 

dr"^ 


+ 1 - 


2M 


-1 


Here we define 


and 


2M di^ (2M ^ l{l + l) 

r2 fj.j. 


2M 
r ' 


-1 


2M\ 

1 - Ip = 0. 


iP 


H ■ 


Then Eq. (Cl) transforms into 

C(1 ~ ~ ~ (^ + l)^w = 0. 

This is the hypergeometric equation. Regularity at r = oo requires 

1 


1p(X-r 


I ’ 


in the limit r —>■ oo. When we require this condition, the unique solution can be written as 


{2M 




/ + !,/ + !, 2 /+ 2 ;- ] , 


where F is the hypergeometric function and q is an arbitrary constant. 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 


APPENDIX D: ASYMPTOTIC EXPANSIONS 


In Eq. (^), if oj ^ 0, r = oo is an irregular singular poin t. A general solution is expressed as a linear combination 
of the purely ingoing and outgoing waves, as seen in Eq. ( |3.lD . The asymptotic expansions for \p^ and \p~ in the 
asymptotic region \uj\r ^ 1 and r ^ M are 


(Dl) 

oo 

V;-(r*)=e-*“’'•^&;r-^ (D2) 

where the coefficients {b^} are given by the following recursive relations: 

(D3) 

^ ^ ^4o;M]6±, (D4) 

[{{v - - l{l + 1) ± iiujM{v - l)}6^_i 

-2M{(!/-1)(2i/-3)-1(/ + 1)}6^_2 (D5) 

+4M^(z/ — 2 )^&y 3 ] for v >S. 
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